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Abstract. Cosmic infrared background (CIB) contains emission from epochs inaccessible to current telescopic studies, such 
as the era of the first stars. We discuss theoretical expectations for the CIB contributions from the early population of massive 
stars. We then present the latest results from the ongoing project by our team @, Q3, Qj], [T^l to measure/constrain CIB 
fluctuations from early epochs using deep Spitzer data. The results show the existence of significant CIB fluctuations at the 
IRAC wavelengths (3.6 to 8 flm) which remain after removing galaxies down to very faint levels. These fluctuations must 
arise from populations that have a significant clustering component, but only low levels of the shot noise. Furthermore, there 
are no correlations between the source-subtracted IRAC maps and the corresponding fields observed with the HST ACS at 
optical wavelengths. Taken together, these data imply that 1) the sources producing the CIB fluctuations are individually 
faint with S v < a few njy at 3.6 and 4.5 flm; 2) are located within the first 0.7Gyr (unless these fluctuations can somehow 
be produced by - so far unobserved - local galaxies of extremely low luminosity and with the unusual for local populations 
clustering pattern), 3) they produce contribution to the net CIB flux of at least 1-2 nW/m 2 /sr at 3.6 and 4.5 flm and must have 
mass-to-light ratio significantly below the present-day populations, and 4) they have angular density of ~ a few per arcsec 2 
and are in the confusion of the present day instruments, but can be individually observable with JWST. 
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Introduction 

Cosmic infrared background (CIB) contains emissions 
produced during the entire history of the Universe, in- 
cluding those arising from the epochs of the first stars 
(see fl5[] for review). Isolating the part of the CIB com- 
ing mostly from the first stars epochs would provide 
important information about the emissions, evolution 
and contents of this era. This paper reviews the recent 
measurements done by our team in deep Spitzer IRAC 
data (Kashlinsky, Arendt, Mather & Moseley - herafter 
KAMM, JKHluKTl) and their interpretation. 

The paper is structured as follows: it starts with a 
brief discussion of the theoretical motivation for our 
measurements, then presents a summary of the most 
recent results at the IRAC Channels 1 (3.6 flm) and 2 (4.5 
jum). Our most recent analysis lfl2ll showed that in the 
fields observed at both IRAC and optical ACS HSTbands 
there are no correlations between the source-subtracted 
IRAC and the ACS source catalog maps. The absence 
of these correlations would place the sources producing 
these fluctuations within the first 0.7 Gyr unless they 
originate in the more local, but extremely low luminosity 
galaxies that somehow escaped ACS detection. Finally, 
we discuss the data requirements on the nature of the 
populations producing these fluctuations: they have to 
originate from very faint sources, located within the first 
Gyr of the Universe's evolution. These require them to 
have mass-to-light ratio significantly below the present- 



day stellar populations and their number density is such 
that they are well within the confusion of the present 
instruments. 



Theoretical expectations 

An early generation of massive stars may have pro- 
duced a non-negligible contribution to the CIB at the 
near-IR wavelengths. The net amount of that contribution 
can account for the entire claimed CIB excess (see re- 
view J5|] for summary of measurements up to 2005) over 
the contribution due to the observed galaxy counts be- 
tween 1 and 4 jim if only ~ 2 — 4% of the baryons went 
through these early stars 16J. The energy spectrum of the 
resultant CIB would be cut off below the Lyman limit at 
the redshift of these stars (~ 1 flm) and the fraction of the 
baryons locked in them would decrease linearly with the 
observed CIB excess flux. The latter would have to be 
below ~ 0.2 — 0.5% if the NIR CIB excess is as low as 
the NICMOS measurements 1 14| suggest. 

The left panel of Fig. Q] shows a template of early 
star emission at z = 10 from lfl3ll . The IRAC filters are 
marked with shaded areas; the figure shows that, al- 
though the IRAC bands are not centered at the peak of 
the emission produced by Ly-a, their wavelengths cover 
the region containing much of the cumulative emission. 
This emission is likely to be reprocessed by the surround- 
ing gas resulting in the model-dependent uncertainty at 
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FIGURE 1. Left: SED from Pop III system is shown for 
z=10. The two solid lines are drawn from [13] for the case when 
re-processing of the radiation takes place in the gas inside the 
nebula and outside in the IGM. Thin dashed-dotted line shows 
the intrinsic SED prior to the emission re-processing by gas. 
Right: Lines show the fluctuations expected for concordance 
ACDM model at z=5 to 20 for the parameters described in the 
text. Shaded area shows the angular scales probed by KAMM. 



least as large as shown in the figure. Note that the color 
of the emission at IRAC channels, while approximately 
flat, can be either slowly or slowly decreasing with wave- 
length. 

The first stars should also have left a distinct imprint 
in the CIB anisotropics, measuring which could give im- 
portant information about that era and its constituents 
md]. There are intuitive reasons why the anisotropics 
can be significant: 1) first stars, if massive, emitted a fac- 
tor ~ 10 5 more light per unit mass than the present-day 
stellar populations, 2) their relative fluctuations would be 
larger because they span a relatively narrow time-span in 
the evolution of the Universe, and 3) they formed (pre- 
sumably) at the high peaks of the underlying density 
field, which would amplify their clustering properties. 

The CIB fluctuations, 8F, on scale 2n/q are related 
to the 2-D power spectrum, P2{q), of the diffuse light 
via 8F ~ [q 2 P 2 {q)/2n] x l 2 . They are related to the CIB 
produced by them, F, and the underlying 3-D power 
spectrum of the sources' clustering, P}(k), via the Lim- 
ber equation which approximates as 8F ~ FA(qd^ (z)), 



where A 2 = — 5" is the mean square fluctuation in 
the source counts over a cylinder of radius ~ k and 
length cAf assuming the sources span cosmic time At. 
The right panel of Fig. 1 shows the typical relative fluc- 
tuation in the CIB produced by the early sources over 
z from 6 to 20 assuming a (realistic) "toy"-model with 
the concordance ACDM initial conditions where the am- 
plification due to biasing is A ~ 0.3(1 + z) y / Af(Gyr). 
The figure shows the distinct peak at ~ 0.2 — 0.3 degrees 
which, for CDM class of models, corresponds to the hori- 
zon at the matter-radiation equality redshifted to z. The 
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FIGURE 2. Shows decrease of the shot noise contribution 
to the power spectrum vs the iteration number of the KAMM 
source subtraction model. Horizontal lines show the levels 
reached in I lOfl . which are a factor of ~ 2 below those in [3]. 
Triangles correspond to the HDFN region, diamonds to CDFS 
region; open symbols to Epoch 1, filled to Epoch 2. 



shaded area shows the range of angular scales covered 
so far by the KAMM measurements. At these scales one 
can expect fluctuations of no more than 3 — 10%, so the 
CIB fluctuation level of 8F ^0.1 nW/m 2 /sr would re- 
quire CIB of at least F ~ (1 -2) nW/m 2 /sr. 



CIB fluctuations from Spitzer deep images 

The KAMM methodology is summarized in I9I [Toh . 
Briefly, the steps are: 1) The maps are assembled using 
the least calibration method of [2]. It solves for detector 
properties at the time of observations and for frame-to- 
frame variations, and does not remove large-scale struc- 
ture except due to (arbitrary) linear gradients. 2) The as- 
sembled IRAC maps are first clipped of resolved sources 
via an iterative procedure which removes all pixels (with 
a certain mask) with flux above a certain threshold. 3) 
The clipped maps are then further cleaned of the re- 
maining emissions from resolved sources using an it- 
erative variant of the clean procedure 0]. 4) The clip- 
ping is stopped so that enough pixels are left (in practice 
> 70%) for robust computation of the fluctuations spec- 
trum via Fourier transform. 5) When clipping is done to 
deeper levels, and progressively smaller fractions of pix- 
els remain, we compute instead the correlation function, 
which is immune to masking [9]. 6) The uncertainties 
on the resultant power spectra (e.g. due to mask which 
destroys the strict orthogonality of the basis functions) 
are verified via simulations. 7) It is verified that, at the 
iterations used for measuring the final fluctuations, the 
subtracted sources components do not correlate with the 
residual diffuse emissions of the maps. 

As of now we have analyzed results from five deep ex- 
posure regions in various parts of the sky |0, [loll . The 
latest results, shown here in Figs. 12131 are from analysis 
of four independent regions mapped in the course of the 



GOODS project observations. They cover ~ 10' on the 
side and have total exposure time of ~ 23 — 25 hours per 
pixel. The power spectrum of the diffuse emission in the 
source subtracted maps is made up of two components: 
1) small scales are dominated by the shot-noise from the 
remaining sources, and 2) larger scales are dominated by 
the component of CIB fluctuations arising from cluster- 
ing of the contributing sources. As we clean maps of pro- 
gressively fainter sources, the shot-noise level decreases 
as shown in Fig. [2] At the same level of the shot noise 
all five regions exhibit the same large-scale fluctuations 
consistent with their cosmological origin. The figure also 
shows the shot-noise levels we stop at when presenting 
the results; they are limited by the floor imposed by the 
instrument noise levels. Comparing the shot-noise levels 
with those computed from the observed galaxy counts 
tells us that individual galaxies have been eliminated to 
AB magnitudes > 26. 

The upper panels in Fig. [3] show the CIB fluctuations 
from the GOODS data at 3.6 and 4.5 fim. Dotted lines 
show the shot-noise. We demonstrated that these fluctu- 
ations do not arise from the various Solar System and 
Galactic foregrounds and are not produced by the in- 
strument systematics HQxJ. The fluctuations also cor- 
relate well with maps at the longest IRAC channels (5.8 
and 8 fim), although at these channels we do not re- 
cover CIB fluctuations with the same S/N due to the 
larger instrument noise and possible cirrus pollution at 
8 jxm. The clustering component dominates fluctuations 
at > 20 — 30'. The solid line shows the fluctuation spec- 
trum due to the concordance ACDM model at high z as- 
sumed for simplicity to have the same amplitude at 3.6 
and 4.5 jim. The figure shows that the signal is well de- 
scribed by the ACDM model with sources at high-z and 
the color of these sources is approximately flat at the two 
IRAC bands. 



Cross-correlating NIR and optical data 
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FIGURE 3. Top panels: CIB fluctuations from the source- 
subtracted maps from 1 12] at the shot-noise levels of 1 10]. Same 
symbol notation as in Fig. [2] Dotted lines show the shot-noise 
contribution. Solid line shows the slope of sources at high-z 
with the ACDM model spectrum of the same amplitude at 3.6 
and 4.5 flm. Lower panels: CIB fluctuations due to ACS galax- 
ies for the 972 x 972 0.6" pixel field at HDFN-Epoch2 region 
for the ACS z-band. Right and left panels show the spectrum 
for the ACS source maps convolved with the IRAC 3.6 and 4.5 
flm beams. Filled circles correspond to ACS galaxies fainter 
than mo = 21 with the mask defined by the clipping. Filled 
diamonds, triangles and squares correspond to fluctuations pro- 
duced by sources fainter than mo + 2, mo + 4, mo + 6. Open 
diamonds show the fluctuations produced by galaxies fainter 
than /;jab=23 when the clipping mask is not applied; the sym- 
bols show that such galaxies were effectively removed from the 
ACS maps by clipping alone done in KAMM. 



The GOODS fields have also been observed with the 
HST ACS instrument at optical wavelengths from ~ 0.4 
to 0.9 pLva reaching source detection levels fainter than 
28 AB mag. If these fluctuations arise from local pop- 
ulations there should be a strong correlation between 
the source-subtracted IRAC maps and the ACS sources. 
Conversely, there should be no such correlations if the 
KAMM signal arises at at epochs where the Lyman break 
(at rest ~ 0.1/im) gets shifted passed the longest ACS 
z-band at ~ 0.9jtim. To test for this we [12] have con- 
structed synthetic maps, overlapping with the GOODS 
fields, using sources in the ACS B,V,i,z bands from the 
ACS sources catalog [3]. These maps were then con- 
volved with the IRAC 3.6 and 4.5 ftm beams. Finally, 



we applied the clipping mask from the IRAC maps and 
computed the fluctuations spectrum produced by the 
ACS sources and their correlations with the IRAC-based 
maps. 

The lower panels in Fig. [3] show the spectrum of 
fluctuations produced by the ACS z-band sources. It is 
clear from the figure that the ACS sources produce a 
completely different spectrum of the fluctuations than 
found by KAMM. This alone shows that they cannot 
produce the KAMM fluctuations signal shown in the 
upper panels. In fact, their power spectrum is very similar 
to that of the CIB fluctuations found in the deep 2MASS 
data, which arise from galaxies at z > 1 |7[]. Fig. 2e of ref. 
||7|] shows the gradual evolution of the power spectrum 
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FIGURE 4. Left: Correlation coefficient between 
clipped/masked ACS and KAMM data. Large and small 
symbols correspond to the IRAC Ch 1 and Ch 2. Open symbols 
correspond to correlations with the maps of the removed 
sources and filled symbols with the residual KAMM maps. 
Circles, diamonds, triangles, squares correspond to HDFN-E1, 
HDFN-E2, CDFS-E2, CDFS-E2. Right: The various thickness 
solid lines show the dimensionless correlation function 
between the diffuse light in the ACS and KAMM maps for 
B,V, (,z-bands. Dotted line shows the dimensionless correlation 
function of the KAMM maps, Ckamm(6)/o"^ amm , which 
remains positive out to ~ 100" and is better viewed when 
presented in log-log plots as in Fig. SI-4 of [9]. 



spectral slope towards deeper cut 2MASS data which is 
consistent with the evolution of the various slices of the 
ACS sources in the lower panels of Fig. [3] 

In ref. Il2ll we compute both the correlation coeffi- 
cients, ,R = (5F acs 5Fkamm)/<7acs<7kamm, and the full 
cross-correlation matrix, C(9) = (8Facs (x)8Fkamm(x + 
9)) between the ACS and KAMM maps. In this repre- 
sentation, the contributions of any white noise (such as 
shot noise and/or instrument noise) component to C{9) 
drop off very rapidly outside the beam and for the IRAC 
3.6 jim channel contribute negligibly to the correlation 
function at 9 greater than a few arcsec. This is best 
seen from the analog of the correlation coefficient at 
non-zero lag, defined as R(9) =C(0)/<7acs<7kamm- The 
mean square fluctuation on scale 9 is given by the inte- 
gral (8F(9) 2 )=jj / e C(9')G'dG'. Hence, we evaluated 

from R(6) a related quantity M(6)=^ J° R(9')9'd9 / , 
shown in the right panels in Fig. [4] 

There are excellent correlations (shown with open 
symbols) between the ACS source maps and the sources 
removed by KAMM prior to computing the remaining 
CIB fluctuations in Fig. [3] However, practically no cor- 
relations remain with the source-subtracted IRAC maps 
and the correlations are negligible outside the IRAC 
beam. This means that, at most, the remaining ACS 
sources contribute to the shot-noise levels in the resid- 
ual KAMM maps (as discussed by flUl ). but not to the 
large scale correlations in Fig. [3] 



The nature of the sources. 

Any interpretation of the KAMM results must repro- 
duce the following major aspects: • The sources in the 
KAMM data were removed to a certain (faint) flux limit, 
so the CIB fluctuations arise in populations with mag- 
nitudes fainter than the corresponding magnitude limit. 

• These sources must reproduce the excess CIB fluctua- 
tions by KAMM on scales > 0.5', both their amplitude 
and their spatial spectrum. • The populations fainter than 
the above magnitude limit must account not only for the 
correlated part of the CIB, but - equally importantly - 
they must reproduce the (low) shot-noise component of 
the KAMM signal, which dominates the power at <0.5'. 

• The data show the absence of correlations between the 
source-subtracted IRAC maps and the same area mapped 
by ACS. • The CIB fluctuations signal is measured from 
the source-subtracted IRAC maps at all four channels 
(although not with the same S/N at the longest wave- 
lengths) and has color which is approximately flat with 
wavelength across the IRAC bands. 

We briefly discuss the constraints in the above order: 

1) Magnitude limits. The nominal limit above which 
sources have been removed in the KAMM analysis is 
wab — 26 at 3.6 and 4.5 ftm and this by itself implies that 
the detected CIB fluctuations arise from fainter systems. 
At this magnitude one is already at the confusion limit 
for the IRAC beam at 3.6 jim, so fainter galaxies can be 
excised only by removing significantly more pixels. This 
magnitude limit corresponds to 6 x 10 h~ 2 L@ emitted at 
6000 A at z=5. A significant fraction of galaxies were 
thus removed from the data by KAMM even at z > 5 and 
the detected CIB fluctuations must be explained by still 
fainter and more distant systems. 

2) Clustering component The relative CIB fluctua- 
tions produced by these sources are at most 10% (see Fig. 
[T]and [11]). Fig. [3] shows that the clustering strength at 
> 1' requires 8Fcm ~ 0.1 nW/m 2 /sr. This meaning that 
the net CIB from sources contributing to the KAMM sig- 
nal at 3.6 jUm is >l-2 nW/m 2 /sr. 

3) Shot noise constraints. The amplitude of the shot- 
noise power gives a particularly strong indication of the 
epochs of the sources contributing to the KAMM signal 
lllll . This can be seen from the expressions for the shot 
noise fsN = ,f( m )Ftot(> mn m ), where f(m) is the flux 
of the source of magnitude m and F tot (> m) is the net 
CIB flux produced by sources fainter than m yD. Fig. 
[3] shows the shot noise amplitude evaluated at 3.6 and 
4.5 jim is P S n — (2,1) x 10~ n nW 2 /m 4 /sr. Above it was 
shown that the sources contributing to the fluctuations 
must have CIB flux greater than a few nW/m 2 /sr and 
combining this with the values for Psn shown leads to 
these sources having typical magnitudes »jab > 29 — 30 
or individual fluxes below a few nJy. 



4) No correlations with optical ACS data. Our re- 
sults show that the source-subtracted CIB fluctuations in 
deep Spitzer images cannot originate in the optical galax- 
ies detected in the GOODS ACS data. While the B,V, i,z- 
band galaxies are well correlated with those seen and re- 
moved in the 3.6 and 4.5 jim data, they correlate poorly 
with the the residual 3.6 and 4.5 jim background emis- 
sion. These galaxies also exhibit a very different spatial 
power spectrum than the KAMM maps and the ampli- 
tude of their fluctuations is generally also low. Whatever 
sources are responsible for the KAMM fluctuations, they 
are not present in the ACS catalog. Since the ACS galax- 
ies do not contribute to the source-subtracted CIB fluctu- 
ations, the latter must arise at z > 7 as is required by the 
Lyman break at rest ~ 0.1 jtim getting redshifted past the 
ACS z-band of peak wavelength ~ 0.9/xm. This would 
place the sources producing the KAMM signal within 
the first 0.7 Gyr. If the KAMM signal were to originate in 
lower z galaxies which escaped the ACS GOODS source 
catalog because they are below the catalog flux threshold, 
they would have to be extremely low-luminosity systems 
(< 2 x 10 7 /z~ 2 L© at z=l) and these galaxies would also 
have to cluster very differently from their ACS counter- 
parts as Fig.|3]shows. 

5) Color. Color gives additional information on the 
nature of the sources producing these fluctuations. Fig. 
|3] shows that the clustering of the populations producing 
the KAMM signal at 3.6 and 4.5 is well described 
by the concordance ACDM model with sources at high z 
with the color approximately expected from populations 
described by SED of the type in Fig.Q] The signal is not 
as clear at the two longest IRAC bands, but the data at all 
four channels correlate meaning that the same sources 
contribute to all the wavelengths [9]. The color at shorter 
wavelengths can prove an interesting diagnostic lfl5ll . 
but requires additional assumptions on the epochs of the 
emissions (which determine the location of the Lyman 
break) and that the differently processed maps leave the 
same populations. 

Additional information on the nature of the popula- 
tions responsible for these CIB fluctuations, can be ob- 
tained from the fact that the significant amount of flux 
(>l-2 nW/m 2 /sr) required to explain the amplitude of 
the fluctuations must be produced within the short time 
available at these high z (cosmic times <0.5-l Gyr). This 
can be translated into the comoving luminosity density 
associated with these populations, which in turn can be 
translated into the fraction of baryons locked in these ob- 
jects with the additional assumption of their F = M/L 
Hill . The smaller the value of T, the fewer baryons are 
required to explain the CIB fluctuations detected in the 
KAMM studies. It turns out that in order not to exceed 
the baryon fraction observed in stars, the populations 
producing these CIB fluctuations had to have F much 
less than the solar value, typical of the present-day pop- 



ulations Hill . This is consistent with the general expec- 
tations of the first stars being very massive. 

For a typical Population III system with an SED such 
as shown in Fig.Q] one expects F ~ 10~ 3 — lO~ 2 r at the 
IRAC bands. Because the (low) shot-noise values at 3.6 
and 4.5 fim imply that the individual sources have flux 
below a few to 10 nJy, the bulk of these populations had 
to have only at most a few times 10 5 M© of stellar mate- 
rial if they are to explain the KAMM fluctuations Hill . 
Such sources would be below the detection threshold of 
the current high-z Lyman-break searches. 

Finally, in order to detect the faint sources responsi- 
ble for the CIB fluctuations with fluxes below a few nJy 
embedded in the underlying sea of emissions, their in- 
dividual flux must exceed the confusion limit 111 111 . If 
such sources were to contribute to the CIB required by 
KAMM data, at 3.6 and 4.5 /im they had to have the 
average surface density of h ~ -Fcib/^sn ~ ^ arcsec -2 . 
In order to avoid the confusion limit and resolve these 
sources individually at, say, 5-sigma level (a = 5) one 
would need a beam of the area (Obeam < 0!~ 2 /n ~5x 
10~ 3 arcsec 2 or of circular radius below ^0.04 arcsec. 
This is clearly not in the realm of the currently operated 
instruments, but the JWST could be able to resolve these 
objects given its sensitivity and resolution. 
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